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ABSTRACT 

We consider here the gravitational collapse of a spherically symmetric inhomoge- 
neous dust cloud described by the Tolman-Bondi models. By studying a general 
class of these models, we find that the end state of the collapse is either a black 
hole or a naked singularity, depending on the parameters of the initial density dis- 
tribution, which are p c , the initial central density of the massive body, and Rq, the 
initial boundary. The collapse ends in a black hole if the dimensionless quantity (3 
constructed out of this initial data is greater than 0.0113, and it ends in a naked 
singularity if f3 is less than this number. A simple interpretation of this result can 
be given in terms of the strength of the gravitational potential at the starting epoch 
of the collapse. 
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The gravitational collapse of a sufficiently massive homogeneous dust ball leads to 
the formation of a black hole, as was indicated by the work of Oppenheimer and 
Snyder [1] , and subsequent investigations. Such a black hole covers the infinite den- 
sity singularity forming at the center of the cloud, which is the end state of collapse. 
However, the final fate of collapse of an inhomogeneous distribution of matter with 
a general equation of state is largely unknown. An important subcase which could 
be examined in this context is the spherically symmetric but inhomogeneous distri- 
bution of dust which collapses under the force of gravity. The general solutions to 
Einstein equations for this case have been given by the Tolman-Bondi spacetimes 
[2] , and it has been demonstrated that naked singularities can occur as the end state 
of such a collapse [3, 4, 5, 6]. In particular, it was pointed out in [5, 6] that the 
collapse ends in a naked singularity or a black hole depending on whether or not a 
certain algebraic equation involving the metric functions and their derivatives has 
positive real roots. 

In order to be able to ascertain the astrophysical implications of such a result, 
it is necessary to translate the condition for the existence of positive real roots into 
the actual constraints on the initial density distribution in the cloud. 

We investigate this issue in the present paper for a class of models of the 
Tolman-Bondi spacetimes. One would expect in this case that the degree of inho- 
mogeneity in the matter distribution plays a role in determining the final fate of the 
collapse. This is demonstrated here by working out the explicit conditions for the 
collapse to the endstate which is either a black hole or a naked singularity depending 
on the initial conditions chosen. It turns out that for the class we are considering, 
these outcomes are characterized in terms of the existence of real positive roots of a 
quartic equation which we shall determine here. This enables us to relate the black 
hole or naked singularity configuration as the end state of the gravitational collapse 
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in terms of the initial density distribution p c and radius Ro of the massive body. 
Finally, we demonstrate our procedure by an explicit example of a typical initial 
density profile from which the collapse develops. 

We use the comoving coordinates (i, r, 6>, 4>) to describe the spherically symmet- 
ric collapse of an inhomogeneous dust cloud. The coordinate r has non-negative 
values and labels the spherical shells of dust and t is the proper time along the 
world lines of particles given by r = constant. The collapse of spherical inhomoge- 
neous dust is described by the Tolman-Bondi metric in comoving coordinates (i.e. 

ds 2 = -dt 2 + ~^jdr 2 + R 2 (t, r) (d6 2 + sin 2 9dcj) 2 ). (1) 

The energy momentum tensor is T* J = e8\8l, where e is the energy density. From 
the Einstein equations it follows that 

' = = m;< « 2 = s + / - (2) 

(We have set 8nG/c 4 = 1). Here the dot and prime denote partial derivatives with 
respect to t and r respectively. The quantity F(r) arises as a free function from 
the integration of the Einstein equations and can be interpreted physically as the 
total mass of the collapsing cloud within a coordinate radius r. Thus we take F > 
for the mass function F(r). The quantity R(ti, r±) denotes the physical radius of a 
shell of collapsing matter at a coordinate radius T\ and on the time slice t = t±. The 
quantities F and / are arbitrary functions of r. In further discussion, we restrict 
to the class of solutions f(r) = 0, which are the marginally bound Tolman-Bondi 
models. Similar considerations could, however, be developed for the models with 
/ > or / < 0. As we are concerned with the collapsing cloud, we take R(t, r) < 0. 

The epoch R = denotes a physical singularity where the spherical shells of 
matter collapse to zero radius and where the density blows up to infinity. The time 



t = t (r) corresponds to the value R = where the area of the shell of matter at a 
constant value of coordinate r vanishes. This singularity curve t = to(r) corresponds 
to the time when the matter shells meet the physical singularity. This specifies the 
range of coordinates for the metric (1) as 

< r < oo, — oo < t < to(r), 

whereas 9 and cf) have the usual coordinate range. In the case of a finite cloud of 
dust, there will be a cut off at r = r&, where the metric is matched smoothly with 
a Schwarzschild exterior. 

The Tolman-Bondi models admit a freedom of scaling in the following sense. 
One could arbitrarily relabel the dust shells given by r =const. on a given t =const. 
epoch, by letting r — > g(r). Thus, at any constant time surface, say at t = to, 
R(r, to) can be chosen to be an arbitrary function of r. This arbitrariness reflects 
essentially the freedom in the choice of units. For convenience of calculation, we 
make the choice of scaling at t = as given by _R(r, 0) = r. With this scaling, the 
R equation in (2), (with / = 0), can be integrated to get 

i? 3 / 2 (r,t)=r 3 / 2 -^v^)t, (3) 

and the energy equation becomes 

^ = ff 2G(r)W, 2G'(r)V (4) 
\ l 3 H(r) J \J 3H'(r)) 

where G(r) = r 3 / 2 , G'{r) = (3/2)r 1 / 2 , and H{r) = y/F(r). We now write F(r) = 
r\(r) and assume A(0) = Ao 7^ and finite, which is the class of models considered 
in Ref. 5. This means that near the origin, F(r) goes as r in the present scaling, 
and the density at the center behaves with time as e(0, t) = 4/3t 2 . This is a general 
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class of models which includes all self-similar solutions as well as a wide range of 
non-self-similar models, which we find quite adequate for the purpose of present 
investigation. The central density becomes singular at t = 0, and the singularity is 
interpreted as having arisen from the evolution of dust collapse which had a finite 
density distribution in the past on an earlier non-singular initial epoch. 

To check if the singularity could be naked, what needs to be examined is the 
possibility that future directed null geodesies could come out of the singularity at 
t = 0, r = 0. The equations of outgoing radial null geodesies in the space-time, 
with k as affine parameter, can be written as, 

d ^- + R'K r K t = (5) 

where K l = dt/dk and K r = dr/dk are tangents to the outgoing null geodesies. 
The partial derivatives R' and R' which occur in (5) and (6) can be worked out 
from Eqn. (3) and are most suitably written as follows, 



R' = T]P — 

A 



2rP 2 



1-r] t 
VA r 

VA r 



R (7) 

(8) 



where we have put 



rF' 

R(r, t) = rP(r, t), r\ = ?](r) = ——. (9) 

F 



The functions r](r) and P(r) have been introduced since they have a well-defined 
limit in the approach to the singularity. 

If the outgoing null geodesies terminate in the past with a definite tangent at 
the singularity (in which case the singularity would be naked), then using (6) and 
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1' Hospital rule we get 

t dt 

X Q = lim - = lim — = lim R' (10) 

i^0,r^0 r t^0,r^0 dr t=0,r=0 

where X = t/r is a new variable. The positive function P(r,t) = P{X,r) is then 
given using (3) by 

2 2P 3 / 2 
X 3~7X ~ ~ 3^ 

We define Q = Q(X) = P(X,0). If future directed null geodesies come out of 
the singularity at t = 0, r = 0, meeting the singularity in the past with a definite 
tangent X = X as given above, then it follows from (11) that such a value X 
must satisfy 



3VA0 

Further, we note, by using the definition F = r\(r), that as r — > 0, r\ — > 1. 
Also, from (2) it follows that for / = 0, UmR = —y^Xo/Q. Using these results in 
the expression (7) for R' implies that the condition (10) is simplified to the equation 

V{X Q ) = (12) 

where 

V(X) = Q + xJ^-X. 

In order to be the past end-point of outgoing null geodesies, at least one real 
positive value of X must satisfy (12) (of course, subject to the constraint implied by 
(11) as stated above, i.e. X < 2/3v^Ao). In general, it is seen [5] that if the equation 
V(Xq) = has a real positive root, the singularity would be naked. Whenever this 
is not realized, the evolution will lead to a black hole. Such a singularity could 
be either locally or globally naked depending on the global features of the function 
A(r). 
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We should clarify the sense in which the terms naked singularity and black 
hole are used. When there are no positive real roots to equation (12), the central 
singularity is not naked, because it follows that there are no outgoing future directed 
null geodesies from the singularity in that case. Further, it is known that the shell 
focusing singularity R = for r > is always covered (for a proof see, e.g. [4], 
[7]). Hence in the absence of positive real roots the collapse will always lead to a 
black hole. On the other hand, if there are positive real roots, it follows that the 
singularity is at least locally naked, though for brevity we have simply called it a 
naked singularity throughout the paper. Such a locally naked singularity would be 
globally naked as well when the outgoing trajectories could reach arbitrarily large 
values of r (i.e. the signals reach far away observers). Otherwise there would still 
be a black hole when these trajectories fall back to the singularity without coming 
out of the horizon. This is a violation of weak censorship only. The occurrence of 
either of these situations will depend actually on the nature of the function A(r). 
The conditions under which this locally naked singularity could be globally naked as 
well have been discussed, for instance, in [6], and we do not go into them here. The 
occurence of positive real roots implies the violation of strong cosmic censorship, 
though not necessarily of weak cosmic censorship. In other words, black hole and 
locally naked singularity are not mutually exclusive alternatives. It can also be 
shown that whenever there is a positive real root to (12), a family of outgoing null 
geodesies always terminates at the singularity in the past [6]. 

We now examine the condition for the occurrence of a naked singularity in 
some detail. Using (11), the condition V(X ) = can be written as 




= 



(13) 




(Recall that F(r) and hence Aq are 
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positive). Using standard results it can be shown that this quartic equation has 
positive real roots if and only if a > a± or a < a2, where 

e*i = — + 5^- 17.3269 
o 

and 

26 r- o 

a 2 = 5V3 « 6.4126 x 10" 3 . 

3 

To derive this condition, first note from (13) that if it has a real root, it must be 
positive as negative values of Y do not solve this equation. Writing the general 
quartic as ax 4 + Abx 3 + 6cx 2 + Adx + e = one defines H = ac — b 2 , 1 = ae — Abd + 
3c 2 , J = ace + 2bcd-ad 2 -eb 2 -c 3 and A = J 3 -27 J 2 . If A < the quartic has two 
real and two imaginary roots. If A > 0, all roots are imaginary unless H < and 
(o 2 / — 12H 2 ) < 0, in which case they are all real. The application of this procedure 
to the quartic in (13) leads to the condition on a given above. 

It follows, however, from (11) as stated earlier that along any such outgoing null 
geodesies we must have yf\^X = Y < 2/3. Then (13) implies that the larger range 
of a for the existence of roots, i.e. a > 17.3269 is ruled out in the sense that no 
outgoing trajectories can meet the singularity with this larger value of tangent Xq. 
This is seen from (13) by writing a as a function of Y, which tells that if a > cki, 
then Y > 2. It thus follows that a naked singularity arises if and only if a < a±, 
or equivalently, if and only if Ao < 0.1809. Whenever the limiting value Ao does 
not satisfy this constraint the gravitational collapse of the dust cloud must end in 
a black hole. The physical interpretation for the quantity Ao can be obtained from 
Eqn.(4) for the time evolution of the energy density. If the collapse starts at a time 
—to < 0, and p c is the initial energy density at the center, then p c = 4/3£q. If p' c be 
the initial density gradient at the center, then we find from (4) that Ao = 16pl/3p' 2 . 

8 



Putting in the units gives that 



A = 



87rG 16p 3 c 
c 4 3p? 



(14) 



Defining (3 = Aq/16, we find that the black hole arises whenever 



0.0113 < (3 



(15) 



and the naked singularity results for the values given by (3 < 0.0113. The occurrence 
of one or the other is governed by conditions on a combination of the initial central 
density and the initial density gradient at the center. The cosmic censorship hy- 
pothesis of Penrose [8] could, in the present context, be translated to the conjecture 
that values of (3 smaller than 0.0113 do not occur in realistic collapse. 

One needs to calculate now the value of the parameter (3 or Ao, given the initial 
density profile for the collapsing massive star. Given the initial central density, one 
can evaluate the initial density gradient dp/dr = p' c at the center as follows. Firstly, 
note that using (4) one can write the expression for dp/dr and it is seen that in the 
limit of r — > this always goes to a finite quantity which is proportional to 1/ \/Ao. 
Now, given the initial data in the form of the density distribution p(R) for the body 
at an initial non-singular epoch of time, in terms of the physical radius R, one can 
integrate using (2) to get the mass function F(R). Then (3) provides a functional 
relationship r(R) which can be inverted (in principle) to express R in terms of r. 
One could then write the mass profile explicitly as F(r) and A is evaluated as the 
limit of F(r)/r as r — > 0. 

These results can be highlighted with an explicit example. To check whether the 
black hole or naked singularity could occur, one will start by examining the initial 
density profile or mass profile for the collapsing body and will go on to calculate A 
as pointed out above. Consider, for instance, the following initial density profile, 
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given at a non-singular epoch t = —t , as a function of the physical radius R: 



e(R, -t Q ) = p{R) 



p c (l- R s /R 3 ), R<Ro 
0, R > R 



(16) 



p c is the central density and Rq the initial boundary of the object. Using (2), it 
follows from the equation for e(r, t) at a constant time epoch that F(r) can be 
written as a function of R as 



(17) 



The relation between R at the initial epoch — to and r can be found from Eqn. (3) 
using the above form of F(R). We get 



r 3/2 = R 3/2 



1/2' 



(18) 



We are interested in finding the form of F(r) near the center, so in (18) we take 
R <C Rq and find after a binomial expansion that R 3 = (16) 1 ' 3 Rlr near the center. 
Using this in (17) shows that near r = 0, the form of F(r) is 



1 



(16) 



1/3 



Rr 



(19) 



Hence, Ao, which is the limiting value of A (r) = F(r)/r, is given by Ao = (16) 1 / 3 p c i?o/3- 
We see that Ao is defined in terms of parameters of the original density profile given 
by Eqn. (16), as expected. After putting in the units, and assuming p/c 2 to be the 
mass density of the collapsing object, and using the condition (15) on /?, we find 
that for the present density profile, given by (16), a black hole will form whenever 
p c Ro/c 2 lies in the range 



1.16 x 10 26 gms/cm < p c R%/c 2 
10 



(20) 



For a comparison, we note that for a neutron star with a central density of 10 15 gms / cm 3 
and radius 10 6 cm, p c R\/c 2 = 10 27 gms /cm. 

We can calculate here the initial density gradient dp/dr = p' c at the center, 
using the form (16) of the density profile, and the relation between r and R near 
r = 0. We get p' c = — (1Q) 1 / 3 p c / Rq which coincides with (14) once we express Ro in 
terms of Ao and p c . Note that in this example, while the physical density gradient 
dp/dR is zero at the center, the coordinate gradient dp/dr at the center is non-zero, 
because the derivative dr/dR at the fixed initial epoch goes to zero at r = 0. In 
fact, we can draw some general conclusions about the density profiles considered in 
this paper. As indicated earlier, p' c goes as l/\/Ao, and is hence non-zero. Using 
(3), we can write, near r = 0, 

dp dp/dR ^Xot^dp/dR 

dr~ = dr/dR = 4R 2 ( > 

Since p' c is non-zero and finite, it follows that dp/dR should behave as R 2 near the 
center. In other words, dp/dR and d 2 p/dR 2 should be zero at r = 0, and d 3 p/dR 3 
should be non-zero. 

Let us define ps to be the value of d 3 p/dR 3 at the origin. It is then straight- 
forward to check for the density profile in (16), as well as for the general density 
profiles considered in this paper, that (3 is proportional to p^ 3 /p^ 3 - Hence it is 
apparent that suitable density gradients, (as measured by ps), will lead, through 
the dynamical time evolution of the collapse, to a naked singularity. On the other 
hand, when p^ is strictly zero, the model is like the Oppenheimer-Snyder model 
for which the singularity is covered. We note, however, that it need not be the 
same as the Oppenheimer-Snyder model, if the higher order derivatives are non- 
zero. Thus we have the situation that while homogeneous collapse leads to a black 
hole, a suitable amount of inhomogeneity, as represented by a non-zero ^3, leads to 
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a naked singularity. One should note that departure from the homogeneous models 
by way of introducing arbitrarily small inhomogeneity does not necessarily make 
the singularity naked, as pointed out by our calculations here. This agrees with the 
earlier results such as those in [3,5,6]. 

Our analysis so far has been exact, within the framework of Tolman-Bondi 
models. While these models may not be sufficiently complex as to describe the 
final fate of real gravitational collapse, they can be useful for the purpose of ob- 
taining insights. With this in mind, we consider the following Newtonian order of 
magnitude estimates. Even though p' c is the central density gradient, defined using 
the coordinate radius r, we approximate it to be the mean initial density gradient 
across the collapsing star. If for an object like a neutron star, we assume the lin- 
ear extent to be L ps 10 6 cm, and assume p c to be (c 2 x mass density), and take 
p c /c 2 10 15 gms / cm 3 , and p' c as p c /L, we get (3 = 0.62. We note that the model 
yields a value of (3 in the range given in Eqn. (15), for which the collapse will lead 
to a black hole, if no force other than gravity is acting. We should emphasize that 
while the condition for the occurrence of the naked singularity is exact when stated 
in terms of the quantity /?, the argument given above is only an approximation and 
is included here as it seems to suggest that the bounds on (3 given here may have 
relevance to real collapse. 

The dimensionless quantity (3 is essentially 2GM / y/Lc 2 , where M = 47rL 3 p m /3 
is approximately the mass of the collapsing body having an initial linear extent 
L and an initial mean density p m , and 7 = p c j ' p m . Looked at in this way, j3/^f is 
nothing but an estimate for the gravitational potential of the body, expressed in the 
basic potential units c 2 /G. It then appears natural that (3 should decide whether 
the collapse ends in a black hole or a naked singularity. Once (3 becomes sufficiently 
large, the collapse appears to proceed to a black hole. From an astrophysical view- 
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point, it is of interest to investigate whether values of (3 large enough to lead to 
a black hole can always be realized in realistic collapse. One also notes that (3 is 
inversely proportional to the degree of inhomogeneity at the center. (This follows 
from (14)). The more inhomogeneous the system, the smaller is f3. 

This model can be generalized in many ways, for instance by considering the 
most general class of functions F(r) and f(r). Also, it has been shown [9] that the 
pattern of a transition from the black-hole configuration to the naked singularity 
configuration persists in models with more general equations of state. It is desir- 
able to cast results for these models in terms of constraints on the initial density 
distribution. Work is in progress to make this connection rigorous. 
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